Shot noise for resonant Cooper pair tunneling 
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We study intrinsic noise of current in a superconducting single-electron transistor, taking into 
account both coherence effects and Coulomb interaction near a Cooper-pair resonance. Due to 
this interplay, the statistics of tunneling events deviates from the Poisson distribution and, more 
important, it shows even-odd asymmetry in the transmitted charge. The zero-frequency noise is 
suppressed significantly when the quasiparticle tunneling rates are comparable to the coherent os- 
cillation frequency of Cooper pairs. 

PACS Numbers: 72.70.-fm, 73.23Hk, 74.50-hr 



Electron tunneling events across a small tunnel junc- 
tion are correlated because of the large charging energy. 
These correlations lead to variety of phenomena which 
fall under the rubric of Coulomb blockade effects [Q. 
As an important example, the single-electron transistor 
(SET) has attracted much interest due to its ultimate 
sensitivity to electric charges If the junctions are 
superconducting, an additional effect, the coherent tun- 
neling of Cooper pairs, comes to play and leads to much 
richer current-voltage characteristics . 

Further understanding of the properties of electron 
transport (related, e.g., to coherence, electron-electron 
interaction and carrier statistics) comes from study of 
current fluctuations In single-electron devices, the 

roles of Coulomb blockade on noise have been discussed 
by many authors |]7|-|ll[ . Moreover the importance of co- 
herence, leading to an enhancement of shot noise in su- 
perconducting quantum point contacts, was pointed out 
in Ref. [|l2|. Up to our knowledge, however, the com- 
bined effect of coherence and Coulomb blockade in su- 
perconducting double tunnel junction systems has not 
been addressed. Additional interest in studying noise in 
single-electron devices comes from their use in quantum 
measurements p^ , p^ a nd as entanglement detectors in 
solid state systems |15[. 

In this Letter, we discuss the statistics of tunneling 
events and the shot noise in superconducting SET near 
a resonance for Cooper pair tunneling. The interplay 
between coherence and interaction, explored by sweep- 
ing the device through the resonance, leads to a num- 
ber of interesting results, i) At a Cooper pair resonance 
the statistics of tunneling events is non-Poissonian and 
it shows an even-odd asymmetry, ii) The shot noise 
suppression depends strongly on the ratio between the 
Josephson coupling and the quasi-particle tunneling rate 
(the effect is more pronounced close to the resonance), 
iii) The frequency-dependent noise has a resonance peak 



at a frequency corresponding to the coherent oscillation 
of Cooper pairs. 

The superconducting SET [see Fig. |l] (a)] is a system of 
two small tunnel junctions in series with a small central 
electrode. The device operates in the regime in which the 
charging energy Ec = e^/2Cs (Cs is the total capaci- 
tance of the island) is much larger than the Josephson 
coupling energy Ej as well as the thermal energy fc^T. 
The largest energy scale is the superconducting gap A 
(assumed equal in both the electrodes and the island). 
By adjusting the bias and gate voltages, V and Vg, one 
can put either the right or left junction at resonance for 
Cooper-pair tunneling We consider the case of 

resonance across the left junction. 

The effective Hamiltonian jl^] is given hy H = Hq + 



\- Ht with y 

Ho^Ecin- 



eVnji — Ej cos 2(j)L 



(1) 



Here n = riL+nn is the number of excess electrons on the 
central island, til {nRj is the number of electrons that 
have passed across the left (right) junction into the cen- 
tral electrode, eno = CbY + CgVg is the offset charge on 
the central island, and 2(/)l is the superconducting phase 



and 



are canoni- 



difference at the left junction, 
cally conjugated [(j)j,nu] = iSjk- ^'-^^ uoii^io ^^pq 
describe the quasiparticles on the electrodes and their 
tunneling across the junctions, respectively Pjl^]. They 
are given by 



Elqp — 
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(3) 



j—L,R kqa 

where 7^^, (jka) creates (annihilates) a quasiparticle with 
momentum k and energy Ska — \/S,'ka + electrode 
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a, ^fe is the single-particle dispersion, and Tkq is the tun- 
neling amplitude. Each event of quasiparticle tunneling 
into {out of) the island across the junctions leads to the 
transition n n + 1 (n ^ n — 1). The rate is given by 



L/R 



{n) = 



Inilqp(g^(f) 

2e 



(4) 



where = ±i?„,„±i, ^ eV ± Enm±i, Em,n = 

Ec{in — n){m -\- n + 2nQ) ^ and Iqp is related to the quasi- 
particle tunneling current pof . 

We will focus on the bias regime \eV\ ~ 2 A + Ec 
Ej, ksT) where two charge states, for example n = 
and n = 2, are nearly degenerate. Then due to the strong 
Coulomb blockade, it suffices to keep the three charge 
states, n = 0, 1,2, and two tunneling rates, Fi = r]^(l) 



and r2 = r^(2); the other tunneling rates are negli- 
gible. To simplify the notation, we will assume that 
Fi = r2 = r, which is a very good approximation in the 
regime we are interested in. Effectively, one can imagine 
that across the left junction only coherent Cooper-pair 
tunneling occurs, interrupted from time to time by quasi- 
particle tunneling across the right junction, see Fig. |l| 
(b). In the experiment of Ref. 1/Fi = 8ns and 

I/F2 = 6ns for Ec = 2.3Ej = 117^eV and A = 2'SOfieV. 

We need to keep track of the variable nn (or alterna- 
tively ul) as well as n {ul = n — nn). Choosing the 
basis of {|n, riij)}, it can be shown that only diagonal 
elements (with respect to nn) of the reduced density ma- 
trix are involved Pmn(nR]t) = {m,nii\ p{t) InjUfi) . The 
generalized master equation can be written in the 

Lindblad form {h = 1): 



n=1.2 



(5) 



where Ln is a Lindblad operator corresponding to the 
quantum jump n — > n — 1, L„ = \n — 1) {n\. The first 
term describes a purely phase-coherent dynamics, while 
the second is responsible for the dephasing and relaxation 
due to the quasiparticle tunneling. 

Counting Statistics. We first investigate the statis- 
tical distribution of the number of electrons that have 
tunneled across the right junction. It has been ob- 
tained p^ , first by defining the characteristic matrix 
G{9, t) — X]„H e^^^"''^p{nR, t + 1), which satisfies a mas- 
ter equation similar to Eq. ^ with the initial condition 
G{0,O) = T,nnPi'^R^i °°)- The probability P(iV,T) 
that N electrons have tunneled during the period r in 
the stationary state is then given by 

d0 



P{N,t) 



27r 



TrG(6',T) 



(6) 



When the dephasing is strong (either F ^ Ej or 
e = Ec[{2+nQf-nl] > Ej), one can show that (TV < 0) 



P(27V,t) 
^ ' ' \N\\ V 2 



\N\ 



exp 



-^f I pa) 



P(27V-1,t) = (7b) 

where F^ = 2£'jF/(4£^ -I- F^) is the relaxation rate for 
the charge state population in the strong dephasing limit. 
The distribution is Poissonian. However, there is a strong 
even-odd asymmetry. Physically, the charge is trans- 
ferred in pairs (i.e. in units of 2e) rather than one by 
one. 

In the weak dephasing limit (F <C Ej) at resonance 
(e = 0), we find 

, Ad 

P{2N, r) = exp 



3Ft\ 



Tdr 



ji^l^l(T), (8a) 



P(2TV-l,T) = -exp 



3Ft\ 



F\N\{r), 



where 



Fn{r) = 

ZTTl 



dz 



■ sinh 



FtA(z) 



(8b) 



(9) 



with \{z) = VI + 82:. This distribution shows a much 
weaker, but still finite, even-odd asymmetry than the pre- 
vious case [see Eq. (|^)]. In the long-time limit (Fr — > 00), 
P(27V,r) = IPg{N,t) and P(27V- l,r) - |Pg(TV,t) 
where Pg{N,t) is a Gaussian distribution with (N) — 
/r/2e and ((AiV)^) = 20Fr/27. 

In the intermediate case (F ~ Ej), an analytic ex- 
pression for P{N, r) is not available. The numerical re- 
sults are shown in Fig. The distribution function de- 
viates significantly from a Poissonian distribution func- 
tion. Coherent oscillations of the Cooper pairs manifest 
themselves in the even-odd asymmetry of the transmit- 
ted charges: P{N, r) is suppressed (enhanced) for odd 
(even) N compared with the Poissonian distribution. 

Shot Noise. The shot noise spectrum is defined as 



S{lu) 



dT e'^^ {{5I{t + T),6I{t)}) , 



(10) 



where 51{t) = I{t) - (lit)) and {A,B} ^ AB + BA. 
The total current /(<) through the system is related to 
the tunneling currents /^/fl — — ec^tn^/fl across the junc- 
tions by 

Ht) = ^lL{t)~^lRit). (11) 

It is convenient to define the spectral densities of tun- 
neling currents Sij{uj) (i,j — L,R) in an analogous way 
as in Eq.(p^ and rewrite the noise power density in the 
form 
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cl 



(12) 



In the stationary state (/) = (1^) = — {Ir)i so that 
S{ijj) = Sll^^^) = SiiR{u!) in the zero-frequency Hmit. In 
the opposite hmit (w oo), S{uj) — {Cl / C^) S fiii{uj) — 
{Cl/CD 2e (/) [In our case, the left junction is 

(nearly) at resonance for the Cooper pair tunneling and 
hence lim(^_>oo Sll{^) = 0; see also the remarks below 
Eq. (§.] 

In order to calculate the two-time correlators in 
Eq. (|l^) , we follow the procedures based on the quantum 
regression theorem starting from the master equa- 
tion (^). An explicit (but lengthy) expression for S{ijj) 
in terms of Fi 2 and e can be given at an arbitrary finite 
frequency (2^. Here we discuss the zero-frequency shot 
noise. 

At = 0, the noise power density takes a simple form 



whose maximal suppression is by factor 1/2 for the 
symmetric junctions. We emphasize, however, that in the 
latter case, the coherence was not essential. In our case, 
the role of coherence is evident noticing that the dip in 
Fano factor [i.e., 5'(0)/2e/] disappears when moving away 
from the resonant condition as shown in Fig. ^. 

In Fig. ^ we show the typical behavior of the finite- 
frequency noise spectrum in the (a) strong and (b) weak 
dephasing limits. It is interesting to notice that (only) 
in the weak dephasing limit, there is a resonance peak at 
Lj = Ej. Near the maximum and for F ^ Ej, the noise 
behaves like |P2| 



5(0) 
2el 



2r2)^2 



(3^5 F2 + 4£2 



(13) 



In the strong dephasing limit (F ^ Ej), the zero- 
frequency shot noise in Eq. (|l^) is enhanced approxi- 
mately by a factor 2 compared with its classical value, 
2el. This can be understood in terms of Josephson quasi- 
particle (JQP) cycle |^,^. Because of the fast quasiparti- 
cle tunneling across the right junction, each Cooper pair 
that has tunneled into the central island breaks up imme- 
diately into quasiparticles, and quickly tunnels out. The 
charge is therefore transfered in units of 2e (compared 
with e in classical charge transfer) for each JQP cycle 
[see also Eq. (^]. 

In the weak (F ^ Ej) and moderate (F ~ Ej) dephas- 
ing limits the semiclassical JQP picture breaks down. In 
the extreme case (F <C Ej), the quasipaticles do not 
see the left junction and consequently the system can 
be viewed (approximately) as a single-junction system. 
Still, the noise deviates slightly from the Poisson value 
since the channels for tunneling (i.e., n = 1 — > and 
n = 2 ^ 1 with corresponding rates Fi and F2) are 
correlated because of the Cooper pair oscillations and 
Coulomb blockade. The effect is related to the resid- 
ual even-odd asymmetry of the distribution function in 
Eq. (|). 

With moderate dephasing (F ~ Ej), quasiparticle tun- 
neling events across the right junction are strongly af- 
fected by the coherent oscillation of Cooper pairs across 
the left junction. Indeed, this effect gives rise to the 
significant deviation from the Poissonian distribution of 
the tunneling statistics. More remarkably, it also leads 
to the suppression of the shot noise which is maximum 
(by factor 2/5) at resonance (e = 0) for F — V^Ej, see 
Fig. H This is reminiscent of the shot noise suppres- 
sion in (non-superconducting) double-junction systems 



5M 
2eJ 



2Cl 



{uj - Ejf + F2/4 



(14) 



The peak is an effect of coherent quantum transitions 
between the two energy levels tunnel split by Ej. 

The JQP process discussed in this Letter was used in 
a recent experiment ||20| to probe the coherent evolution 
of quantum states in a Cooper pair box. Weak continu- 
ous measurement using quantum point contact [ psf and 
strong measurement using single-electron transistor [T^ ] 
have been proposed. Whereas both schemes are non- 
invasive measurements, the setup discussed here probes 
the charge states on the island directly and invasively. In 
the weak dephasing limit, the resonance peaks in Eq. ( p^ ) 
and in Ref. ]2^ have a similar physical origin, and yet 
the latter has a peculiar upper bound. In the strong 
dephasing limit, the broad peak around zero frequency 
(a)] does not fit to the single Lorentzian shape of 
|13| , which is a manifestation of the random tele- 
graph noise; nevertheless, F^ and l/rmix in Ref. give 
the same time scale describing the relaxation of popula- 
tion density of the charge states on the island. 

In conclusion, we have investigated the combined ef- 
fects of coherence and interaction on the statistics of 
tunneling events and the shot noise in a superconducting 
SET. It has been shown that the number distribution of 
tunneled electrons deviates from the classical Poisson dis- 
tribution and that zero-frequency shot noise is suppressed 
significantly due the coherent oscillation of Cooper pairs 
in the presence of Coulomb blockade. 
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invaluable discussions. We also acknowledge useful 
discussions with Y. Blanter, C. Bruder, G. Falci, Y. 
Makhlin, A. Maassen van den Brink, G. Schon, and J. 
Siewert. 
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FIG. 1. Schematic diagrams of (a) the superconducting 
SET device and (b) the transition processes between relevant 
charge states. 




FIG. 2. Probability distribution function P(N, r) at 
Ft = 4 for F = y/2Ej (solid line). For a comparison, the 
Poissonian distribution is also plotted (dotted line). Notice 
that TV < by definition. 
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FIG. 3. Normalized zero-frequency shot noise for 
e/Ej — 0, 0.25, • ■ • , 5. The dip in the noise is most pro- 
nounced at resonance (e/Ej — 0). 
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FIG. 4. Typical behavior of noise power spectrum S{u)) as 
a function of frequency lo in the (a) strong (ri,2 ^ Ej) and 
(b) weak (ri,2 Ej) quasi-particle tunneling limits. For 
both plots, Cl ~ Cr — Ce/2 were assumed. 



